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■ In a previous paper we introduced examples of Hamiltonian mappings with phase 

space structures resembling circle packings. We now concentrate on one particular 
mapping and present numerical evidence which supports the conjecture that the set 
of circular resonance islands is dense in phase space. 
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I. INTRODUCTION 

Piecewise isometries are iterated mappings which preserve distances on each of a number 
of disjoint regions of space. Such maps and their variants are currently enjoying considerable 
attention in the literature [I], 0, |J || |[ [7|, |8], |9|, [10], [11], fT2| , and arise naturally in the 



context of digital filters \SQ g§ g§ %Q f§ (Eg, [2g, |2IJ and dual billiards |22], |23], |24f. In 
a recent article |25[] three new piecewise isometries were introduced, each corresponding to 
a particular phase space geometry: planar, hyperbolic or spherical. As examples of kicked 
Hamiltonian systems, all of these maps invite physical interpretations. For instance, in the 
presence of damping, the planar map describes oscillations in a water-filled U-tube under 
(Nj ■ the periodic kicking of bubbles [^(J . 

In the present paper we restrict our attention to the mapping defined on a sphere. It has 
a dynamical behaviour which is dominated by an infinite number of stable periodic orbits of 
arbitrarily long period. The resonance islands associated with these orbits are circular discs 
which fill in the phase space in a manner resembling circle packings |27], |28, 29], |30| . 
. ^ | A circle packing is a set with empty interior and whose complement is the union of disjoint 

open circular discs. Circle packings with zero Lebesgue measure are called complete. We 
have previously conjectured f25[ that the phase space structure produced by this map is 
indeed a circle packing, but unlikely to be complete. We now present numerical evidence 
which supports this belief. Our choice to concentrate on the spherical map was due to the 
compactness of its phase space. This proves to be of a great advantage in the numerical 



computations. However all three of the maps described in |25[] are conjectured to produce 
circle packings. 

Numerical studies have led Ashwin et al jl], Q to put forward similar conjectures for the 
sawtooth standard map (see also Goetz for proofs). This mapping might be considered 
an analogue of those described above for the torus. However the fundamental question about 
the 'packed' nature of the phase space remained untested. To this endeavour the notion of 
riddling was developed []TT| [TT], Also of interest in the literature is the concept of a fat 
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fractal M, M [35L M, B71 B8L BS 



Fat fractals are loosely defined to be sets with a fractal 
structure but non-zero measure. The circle packings described above might be considered 
examples of such sets. 

To test whether the set of circular resonances is dense in phase space, we will employ 
a symbolic labeling of the periodic orbits and use special symmetries in the mapping, to 
efficiently calculate a large class of the discs. Our numerical studies indicate that this class 
may in fact itself pack the phase space, and therefore, is the entire set. 

Our paper is organized as follows. In the next section we will give an overview of the 
mapping to be analyzed, stating three main conjectures concerning its behaviour. In section 



III] we present numerical evidence which supports these conjectures. Finally, in section [TV 



we will discuss our results. 



II. THE MAP 



The Hamiltonian under consideration is that of a kicked linear top 

oo 

H(J,t) = ujJ 3 + fi\Ji\ 6 ( f 



n) 



(1) 



where [A, uj G [0, 2ir) are parameters and (J)j = J, = tijkXjPk {i = 1,2,3) are the three 
components of angular momentum for a particle confined to a sphere, normalized such that 
J ■ J = 1. The evolution of J is governed by the equations 



where {■ , ■} are the Poisson brackets. Their solution can be written as the mapping 
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where s n = sgn J™, which takes J from just before a kick to one period later. Here sgnx is 
the signum function with the convention sgnO = 0. 

The mapping is deceptively simple. It first rotates the eastern hemisphere (J\ > 0) 
through an angle fi and the western hemisphere (Ji < 0) through an angle — fi. The great 
circle J\ — remains fixed. Then the entire sphere is rotated about the J3 axis through 
an angle uj. Thus the map rotates every point on the sphere in a piecewise linear fashion 
except on Ji = where its Jacobian is singular. Surprisingly, this leads to a phase space of 
great structural complexity. An example is plotted in figure [y. The eastern hemisphere is 
shown in black, the western in gray. To simplify our study we will restrict ourselves to the 
case fi = uj = ir(\/5 — 1), as shown in the figure. 

Note that if an orbit of our map does not have a point on J\ = 0, then its Lyapunov 
exponents are zero, making it stable. If however, it does have a point on J\ = 0, then 
its Lyapunov exponents are undefined. We overcome this problem by simply defining an 
unstable orbit to be one with a point on J\ = 0, and all others stable. The unstable set is 
defined to be the closure of the set of all images and preimages of the great circle J\ = 0. In 
figure |l] the unstable set is in black and gray. The circular holes are resonances consisting 



FIG. 1: The phase space when \i = u = 7r(\/5 — 1). 



entirely of stable orbits. At their center lies a stable periodic orbit. All other points inside 
the resonances rotate about the central periodic orbit in a linear fashion. The period of 
rotation is an irrational multiple of 7r for almost all values of the parameters /i and u, hence 
their circular shape. If the period of rotation happens to chance upon a rational multiple of 
Ti, then the resonance forms a polygon (see |25|, |40[|). 

We have previously shown |25| that it is possible to label every point of a periodic orbit 
of least period n with a unique sequence {s k = 0, ±l}fc =0 .. n -i- The position of the point is 
given by the solution of 

J = F(s n - 1 )F(s n - 2 )...F(s°)J = RJ (4) 

lying on the unit sphere with sgn J% = s°. The other n — 1 points of the periodic orbit are 
found by substituting the n — 1 cycles of the original sequence into (||). Hence, if {J fc }fc=o..n-i 
is a periodic orbit then the sequence corresponding to its first point is {sgn J\ }k=o..n-i- The 
orbit is unstable if and only if this sequence contains 0. Note that the (unnormalized) 
solution of (^) is simply the axis of rotation of R, namely 



J = (R23 — R32 , R31 — R43 > R12 — R2l)- 



(5) 
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Although every periodic point is uniquely represented by a sequence, not every sequence 
represents a periodic point. Hence we still need to find which sequences are legitimate. This 
may be done by checking each of the 2" different possible sequences for a period-n orbit. 
However, this is computationally expensive for large periods. An alternative approach is to 
exploit symmetries of our mapping. 



If we decompose our map into involutions f41| 

F = I 2 o Ji = I 2 o I x 



(6) 



where 



iiJ = {-J\ , J 2 cos(/i sgn Ji) - J 3 sin(/i sgn Ji) , J 2 sin(/i sgn Ji) + J 3 cos(/x sgn Ji)) (7) 

J 2 J = (— Ji cosco> — J 2 sinu; , — Ji sinu; + J 2 coscj , J 3 ) (8) 

Ji J = ( Ji , - J 2 cos(jU sgn Ji) + J 3 sin(/i sgn Ji) , J 2 sin(/i sgn J x ) + J 3 cos(^ sgn Ji)) (9) 

J 2 J = ( Ji cos uj + J 2 sin uj , Ji sin uj — J 2 cos , J 3 ) (10) 



and 
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(11) 



we obtain the four symmetry lines 



Si 


Jl 
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(12) 


s 2 


Jl 


= — J 2 tan(u;/2) 


(13) 


Si 


J 2 


= J 3 tan(/i/2) sgn Ji 


(14) 


s 2 


J 2 


= Jitan(a;/2) 


(15) 



as the fixed lines (Fix(J) = {x\Ix = x}) of our involutions. Given that IFI = F 1 for all 



the above involutions, it is not difficult to prove the following theorem [33 



Theorem. Let {J fc }fc=o..n-i be a periodic orbit of F with least period n. Suppose 
that the point J° lies on a symmetry line of F (Si, S 2 , Si or S 2 ). Then no other points of 
the orbit lie on the symmetry line if n is odd, or exactly one other, J n//2 ; if n is even. 

This property reduces the task of finding periodic orbits to a one-dimensional search. 
However, if we wish to find all periodic orbits of a given period, then we must first establish 
that every periodic orbit must have a point on at least one of the four given symmetry lines. 
By definition, every unstable orbit has a point on Si, and thus, we have the following fact. 

Fact. Every unstable periodic orbit has a point on Si. 



In [25 we made the following conjecture based on preliminary numerical work. 



Conjecture 1. For almost all u and ji, every stable periodic orbit has a point on 
at least one of S 2 , Si or S 2 . 



The conjecture is not true for all u and /i since, as remarked above, the rotation of 
a resonance may become a rational multiple of 7r when the parameters take on particular 
values. However this occurs on a set of measure zero. 
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The above theorem, and the corresponding symmetric periodic orbits for which it refers, 
have become folklore in the theory of reversible dynamical systems [52]. Proponents of sym- 
metry methods include Birkoff El, De Vogelaere [0] and Greene et al. ||45|| . Unfortunately, 
one generally finds that many other periodic orbits will be present which do not have points 
on a symmetry line. Thus, the mapping we are considering is conjectured to be of a very 
special class for which all periodic orbits are of the symmetric type. 

We now turn our attention to the unstable set. The complicated structure of circular 
resonances with decreasing radii displayed in figure [I] suggests the following conjecture. 



Conjecture 2. The unstable set has empty interior. 



That is, the unstable set is a circle packing of the unit sphere. We also put forward 
a third conjecture. 

Conjecture 3. The unstable set has positive Lebesgue measure. 

Or equivalently, the circle packing is not complete. Numerical evidence which sup- 
ports the above three conjectures will now be presented. 



III. NUMERICAL EVIDENCE 



Suppose we ignore the fact that the stable symmetric periodic orbits may only form a 
subset of the set of all stable periodic orbits. If we find that the set of circular resonances 
associated with these orbits pack the phase space (i.e. the set of resonance discs is dense on 
the surface of the unit sphere), then we are left with only one conclusion: both Conjecture 
1 and 2 must be correct. Furthermore, we are now in a position to test Conjecture 3. This 
is precisely the approach we will take. Consider the following trivial property. 

Suppose J° lies on a symmetry line corresponding to the involution / (JJ° = J°), and 
the iterate J m also lies on this same symmetry line (IF m J° = F m J°). Then, 

F 2m J° = F m F m J° = F m IF m J° = IF- m I 2 F m 3° = 13° = J° (16) 

and hence, J° must be a periodic point with period 2m. To find the orbit's least period we 
construct the unique sequence {sgn J± }k=o..2m-i and locate the shortest word which generates 
the orbit. If J m was the first iterate to hit the symmetry line, then either the period of the 
orbit is odd with period m, or even with period 2m. 

Thus we only need to find where each symmetry line intersects with its image; an effortless 
task to solve numerically by virtue of the simplicity of our mapping. Every image of a 
symmetry line is just a collection of arc segments on the sphere. 

Now suppose that we have just found a particular periodic orbit. If we are to gain some 
insight into the above three conjectures we will need to find the size of its circular resonances. 
Note that the resonance of at least one point of the orbit must touch the great circle Ji = 
(Si). This is because the boundary of each resonance forms part of the unstable set, and 
hence, iterates arbitrarily close to S\. Consequently 

\Ji\ (17) 



a = mm 

k=0..n- 



where n is the period of the orbit, and a is perpendicular distance between the boundary of 
the circular resonance, and the line connecting its center (the periodic point) with the origin. 
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FIG. 2: (a) The radii of the circles for increasing periods of the orbit, (b) The number of circles 
with radius > r. (c) The residual area not covered by the set of circles with radius > r. (d) The 
radius of the largest circle which will fit into the complement of this set. 



All resonances of the same periodic orbit are of equal size. The sections of the symmetry line 
which intersect the orbit's resonances contain no other periodic orbits and may be ignored 
for the remainder of our search. 

For the symmetry lines 5*2, Si and £>2, the above procedure was run for m = 5 x 10 5 
iterations of the mapping. Hence all symmetric periodic orbits with odd periods < 5 x 10 5 
and even periods < 1 X 10 6 were found, producing over 22 billion circles. 

To investigate properties of the unstable set we will define the radius of a circle on the 
unit sphere to be 

r = 2 sin(arcsin(cr)/2) (18) 

which allows us to write the spherical area covered by the resonance as nr 2 . Other definitions 
for the radius such as a, or the angular radius arcsin(o"), are also suitable since all are 
equivalent in the limit a — > 0. 

The proliferation of circles with smaller and smaller radii is displayed in figure |2|(a). Here 
the radius of the resonance versus the period of the orbit is plotted on a logarithmic scale. 
Note that we can only be sure that we have all circles of a particular radius for radii ^ 10~ 5 ' 5 
(the upper dashed line in figure @(a)). 



7 



Next we define the quantity N(r) as the number of circles of radius larger than r. It is 
found that N(r) follows the scaling law 

N(r) = Ar-- d (19) 

(see figure 0(b)). Linear regression on a data set in the range 10~ 4 < r < 10~ 5 ' 5 reveals the 
exponent to be d ~ 1.8. If the unstable set is indeed a circle packing, then the above scaling 
law implies that d will be its circle packing exponent 

d = sup <x:y^rf = oo> = inf < y : r ? < oo > (20) 

where the radii rj are labeled to form a non- increasing sequence: r± > > The circle 

packing exponent is known as the Besicovitch- Taylor index in one dimension | 46|| , and is 



equivalent to the exterior dimension d x and uncertainty exponent a of Grebogi et al [35 



d = d x = 2 — a 37 



To investigate the Lebesgue measure of the unstable set we define the residual area 

R(r)=A7r-7rJ2 r ] ( 21 ) 

3=1 

where j(r) is the largest j such that rj > r. This quantity follows the scaling law 

R(r) = R + Br 2 - d . (22) 

Using linear regression on our data set (see figure 0(c)) we find that Rq/^tt « 0.1169, and 
hence, the set of resonances (for the symmetric periodic orbits) measure less than 89% of 
the sphere. Note that the above scaling law also implies that 2 — d = (3, where (3 is the 



fatness exponent of Farmer |33 



Finally we ask the most critical question: Does the set of resonances for the symmetric 
periodic orbits pack the surface of the sphere? A numerical scheme for answering this 
question might be to consider the set C r of all circular resonances with radius larger than r. 
If we define the quantity sir) to be the radius of the largest circle which will 'fit' into the 
complement of this set, then the unstable set has empty interior if and only if s(r) — > as 
r — > 0. In figure 0(d) we plot s(r) on a logarithmic scale and find that the scaling law 

s ( r ) = Cr p (23) 

is approximately obeyed. The exponent might be thought of as describing the 'completeness' 
of the circle packing, and invites the general definition 

^ lim lo^(r)_ 

r^O log r 

When constructing Apollonian packings, s{r) is the radius of the next circle to be packed, 
and hence p — 1. In our case p ~ 0.6. This estimate, however, might prove to be a poor 
one. Note that s(r) has been calculated only for r < 10 -4 . Ideally we would like to use our 
entire set of circles and increase the graph to r = 10 -5 5 . However the current plot took over 
a year to produce. Our procedure for finding all symmetric circles with radius larger than r 
is extremely efficient. However we are unaware of an efficient algorithm for the calculation 
of s(r). 
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IV. CONCLUSION 



In this paper we have outlined a procedure to efficiently calculate all symmetric periodic 
orbits up to a given period for the spherical map of |5|. When // = u = n(y/E - 1), the 



circular resonance islands corresponding to the stable orbits are found to initially follow a 
scaling law which indicates that they 'pack' the surface of the unit sphere. This would imply 
that all periodic orbits are of the symmetric type. Under this assumption, the unstable set 
has empty interior, but appears not to be of zero Lebesgue measure. The latter result 
agrees with Ashwin et al [|1], EJ for the sawtooth standard map. As a final remark, it might 
prove fruitful to study the index p (EMI) as an indicator for nowhere dense sets. Measures of 
connectedness have already been defined K7 . 
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